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ABSTRACT 



Aims. We define the relationship between the double-mode pulsation of Cepheids and metallicity in a more accurate way, determine 
the empirical metallicities of double-mode Cepheids from homogeneous, high-resolution spectroscopic data, and study of the period- 
ratio - metallicity dependence. 

Methods. The high S/N echelle spectra obtained with the FEROS spectrograph were analyzed using a self-developed IRAF script, 
and the iron abundances were determined by comparing with synthetic spectra assuming LTE. 

Results. Accurate [Fe/H] values of 17 galactic beat Cepheids were determined. All these stars have solar or slightly subsolar metallic- 
ity. Their period ratio (Pi/Po) shows strong correlation with their derived [Fe/H] values. The corresponding period ratio - metallicity 
relation has been evaluated. 

Key words. (Stars: variables:) Cepheids - Stars: atmospheres - Stars: general 
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1. Introduction 

It was realized decades ago, that the period ratio of double-mode 
Cepheids (i.e. classical Cepheids that are pulsating in two modes 
simultaneously, also called beat Cepheids) strongly depends on 
the physical parameters of these stars, namely the mass (M), lu- 
minosity (L), mean effective temperature (T e ff), and metallicity 
(Z). This can be used to give tight constraints on these phys- 
ical parameters via comparison of observed periods and pre- 
diction of theoretical pulsational models on the Petersen dia- 
gram {P\ I P a vs. logP ; Petersen [1973] M organ & Welch [19971 
D'Cruz, Morgan, & Bohm-Vitense 2000). The comparison of 
beat Cepheids in both Magellanic Clouds (LMC and SMC) with 
their galactic counterparts revealed that the key parameter is the 
metallicity. In order to get unambiguous estimates for the mass 
and luminosity of a given Cepheid, its metallicity must be known 
from observations. 

Besides constraining pulsational and evolutionary models, 
metal abundances of Cepheids are useful for studying the abun- 
dance distribution in the Milky Way (Kovtyukh et al. 2005; 
Lemasle et al. I2007l l. as well as in other galaxies (Beaulieu et 
al. [2006). Recent theoretical calculations also indicate the im- 
portance of metallicity for the double-mode pulsation (Buchler 
& Szabo|2007]l. 

Unfortunately, double-mode Cepheids are not as frequent as 
their single-mode counterparts, because radial pulsation with si- 
multaneously excited two periods can only be maintained in a 
very limited region of the instability strip. The number of known 
beat Cepheids in the Milky Way galaxy is still very small, only 



Send offprint requests to: K. Sziladi 

* Based on observations taken with the ESO 2.2-m telescope at La 
Silla, Chile (Proposal 073.D-0072) and the 1.82 m telescope at David 
Dunlap Observatory, Canada. 



23. The situation is better for LMC and SMC, where we have 
better statistical samples of the beat Cepheid population, thanks 
to the long-term microlensing surveys. 

In this paper we present new iron abundances for an almost 
complete sample of known double-mode Cepheids in the Milky 
Way. These data are based on analyzing of high-resolution, high 
signal-to-noise echelle spectra. The observations and the data re- 
duction are described in Sect. 2. In Sect. 3 we give the details of 
the applied method of abundance analysis. The results are pre- 
sented and discussed in Sects. 4 and 5. Finally, we draw our con- 
clusions in Sect. 6. 

2. Observations and data reduction 

The observations were carried out with the fiber-fed high- 
resolution echelle spectrograph FEROS at the ESO/MPG 2.2-m 
telescope at La Silla observatory, Chile (Kaufer et al., 1999). 
FEROS has a continuous wavelength coverage from 3900 to 
9200 A. Two fibers simultaneously recorded the light from 
the object and the nearby sky background. The detector was a 
back-illuminated CCD with 2148x4096 pixels of 15 /mi size. 
Seventeen program stars were selected as observing targets from 
the list of known galactic double-mode Cepheids accessible 
from the southern hemisphere (Table 1). No previous extensive 
investigation of double-mode Cepheids has ever been performed 
at the resolution achieved by FEROS (R = 48000). 

The spectra of the program stars were collected on 3 nights 
between 30 May and 1 June 2004. The observed frames were 
read out in the "normal" mode with readout time ~41.6 s, detec- 
tor gain ~3. 2 e /ADU, and readout noise ~5.1 e~. The exposure 
times were chosen based on primarily the brightness of the tar- 
get star, but avoiding overly long exposures during which the 
stellar parameters may vary significantly due to the pulsation. 



2 



K. Sziladi et al.: Iron abundances of double-mode Cepheids 



Table 1. Basic parameters of program stars. References: A: 
Antipin (fT997l ; BP: Beltrame & Poretti (|2002j l; PP: Pardo & 
Poretti <fl997V , WO: Wils & Otero ( f2004b . 



Star P a (days) P\/Pq Reference 



Y Car 


3.6398 


0.7032 


PP 


EY Car 


2.8760 


0.7079 


WO 


GZ Car 


4.1589 


0.7054 


PP 


TU Cas^ 


2.1393 


0.7097 


PP 


UZCen" 


3.3343 


0.7064 


PP 


BKCen 


3.1739 


0.7004 


PP 


VX Pup 


3.0118 


0.7104 


PP 


EWSct 


5.8235 


0.6986 


PP 


V367 Set 


6.2932 


0.6968 


PP 


V458 Set 


4.8412 


0.6993 


A 


BQ Ser 


4.2710 


0.7052 


PP 


UTrA 


2.5684 


0.7105 


PP 


AP Vel 


3.1278 


0.7033 


PP 


AX Vel 


3.6732 


0.7059 


PP 


GSC 8691-1294 


4.317 


0.7035 


WO 


GSC 8607-0608 


4.089 


0.7017 


wo 


VI 048 Cen 6 


0.9224 


0.8058 


BP 



" data obtained at DDO 
6 lst/2nd overtone pulsator 

Calibration data, including bias, dark, flat-field, and ThArNe 
lamp exposures were obtained during the standard calibration 
routine. Moreover, several daytime solar spectra were recorded 
with the same setup before each night. 

In addition to the program stars observed with FEROS, 
two echelle spectra of the northern beat Cepheid TU Cas were 
also included in the sample. These spectra were taken at the 
David Dunlap Observatory, Canada, with the echelle spectro- 
graph mounted on the 74 inch (188 cm) telescope in 1997. See 
Kiss & Vinko (2000) for details on the instrument and the reduc- 
tion process. 

The reduction of the FEROS data was done in IRAF, inde- 
pendently of the ESO/MIDAS pipeline. We preferred an inter- 
active reduction, in order to check each reduction step carefully 
and avoid any possible misinterpretation of the resulting spectral 
features. 

First, all frames were bias-corrected by subtracting a master 
bias frame, defined as the average of several bias frames corre- 
sponding to the "normal" readout mode (see above). The stel- 
lar frames were then corrected for cosmic rays with the task 
noao.imred.crutil.cosmicrays. Then, an averaged master flat- 
field frame was created in the same way as the master bias frame. 
The echelle orders were identified on this master flat-field frame. 
Since the echelle orders showed complicated double-peak flux 
profiles that had a varying shape across the chip, the automatic 
order- finding algorithms in IRAF could not provide reliable re- 
sults. Thus, we visually identified all spectral orders and manu- 
ally adjusted the position and width of each aperture centered on 
the particular order. This way we could reliably define and ex- 
tract the brightest 35 orders both for the object and the sky fiber, 
respectively. 

After this step, we independently followed two different 
approaches during the reduction. For the first, we applied 
the noao.imred. echelle. dofoe task defined particularly for han- 
dling two-fiber-fed echelle spectra. This complex task performs 
scattered-light subtraction, order extraction based on an order- 
reference frame, lD-flatfielding using the extracted flat-field or- 
ders, and wavelength calibration based on a spectral lamp frame. 



The order-reference frame was the master flat-field frame de- 
scribed above. The extracted flat-field orders were fitted by a 
smooth polynomial, keeping only the pixel-to-pixel variations, 
and these "residual" flat-fields were used in the flat-field divi- 
sion. This kind of flatfielding does not correct the stellar spec- 
tra for the blaze function, which must be determined and re- 
moved subsequently. This was done along with the continuum 
normalization later. For the wavelength calibration, a master 
arc frame was created by median-combining several long- and 
short-exposure arc frames obtained during the calibration pro- 
cedure. The reason was that the short-exposure arc spectra had 
poor signal-to-noise (S/N) for many weak lines, while the long- 
exposure ones had too many badly saturated lines. The combi- 
nation of the arc frames was tuned to have as many weak lines 
with good S/N as possible with the smallest number of saturated 
strong lines. The stability of FEROS enabled us to use this single 
master frame for the wavelength calibration of all object spec- 
tra obtained on the same night. The wavelength calibration was 
performed within dofoe by a two-dimensional polynomial fit- 
ting to the order numbers and the measured pixel positions of 
the arc emission lines. Finally, the resulting object spectra were 
continuum-normalized by fitting a smooth polynomial to the lo- 
cal continuum level, determined iteratively by rejecting those 
pixels that deviated more than 2cr from the local mean level. 

The second approach was based on a self-developed CL- 
script that includes other IRAF routines for reducing and cali- 
brating general echelle spectra. The reduction steps were basi- 
cally the same as described above. Specifically, order extraction 
was done with the noao. imred. echelle. apall task. Instead of fit- 
ting a global illumination pattern for the scattered light, a local 
background flux level was determined and subtracted from each 
order during the extraction. The flat-field correction was done by 
dividing each extracted order of the object frame by the corre- 
sponding extracted flat-field order. This way the blaze function 
was also removed from the object spectra, but the spectrum of 
the flat-field lamp distorted the stellar (and background) fluxes. 
However, since our purpose was to get continuum-normalized 
spectra, this did not cause any problem, because the final con- 
tinuum normalization easily removed the slight flux tilt caused 
by the flat-field spectrum. In fact, we found that this approach 
produces better normalized stellar continuum than the first ap- 
proach described above, because it avoids the numerical prob- 
lems encountered in the simultaneous fitting of the strongly 
varying blaze function and the stellar continuum. As in the 
first method, the wavelength calibration was based on a two- 
dimensional polynomial fit to the extracted orders of the mas- 
ter arc frame, using the task imred. echelle. ecidentify. For the sky 
apertures, the task imred. echelle. ecreidentify was applied based 
on the wavelength solution of the corresponding object aper- 
ture. As in dofoe, this accounts for the slight zero-point shift 
between the wavelength solutions of the corresponding object 
and sky apertures. The removal of the sky spectra from the stel- 
lar spectra was done on the wavelength scale, instead of pixel 
scale, with the task imred.echelle.sarith. The length of each or- 
der was then manually adjusted by trimming the noisy edges at 
the ends. Finally, all object spectra were continuum-normalized 
using noao. imred.echelle. continuum, as described above. 

The final spectra from the two independent reductions were 
compared order-by-order, and very good agreement was found. 
A few reduced spectra were also double-checked with those 
obtained by the imporved FEROS pipeline in MIDAS (Rainer 
2003 ). Very good agreement was found, except for a slight wave- 
length shift, which is due to the fact that the MIDAS pipeline 
contains a barycentric Doppler-correction of the final spectra. 
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Fig. 1. Fe abundances for BQ Ser with optimal atmospheric pa- 
rameters. The parameters were set as r e ff = 6000 K, log g = 2.0, 
Vt = 3.7 km s _1 . Open circles denote Fe I lines, filled squares are 
Fe II lines. Top panel: iron abundance vs. excitation potential; 
bottom panel: abundance vs. EW. 



The precision of our wavelength calibration was tested by com- 
paring the wavelengths of the telluric features in the vicinity of 
Ha with those from other spectra obtained with other instru- 
ments. Perfect agreement, within 1 pixel (~0.05 A), was found. 
We believe that our reduction efforts have resulted in reliable, 
properly wavelength-calibrated, continuum-normalized FEROS 
spectra of all our program stars. Our calibration also supplies 
a reduction approach within the IRAF environment to potential 
users as an alternative to the standard MIDAS pipeline (Sziladi, 
PhD thesis, in preparation). 



3. Abundance analysis 

In order to determine the iron abundance of the program stars, 
first, a list of weak, unblended Fe I and Fe II lines was assembled 
based on the information (wavelength, oscillator strength, and 
excitation potential) found in the extensive line list of Kovtyukh 
& Andrievsky (1999). Each line profile was visually examined 
and blended ones were rejected. This resulted in 77 Fe I and 
18 Fe II lines. Then, the equivalent widths (EWs) of these lines 
were calculated from the observed spectra using the IRAF task 
splot. First, the spectrum was Doppler-corrected (with the task 
dopcor) to bring the lines as close to their laboratory wavelength 
as possible. The input radial velocity was calculated via cross- 
correlating the stellar spectrum with a daytime solar spectrum, 
which has approximately zero radial velocity with respect to the 
observer. The 20th echelle order containing the Mgl 'b' triplet at 
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Fig. 2. The sensitivity of the calculated line abundances to the 
atmospheric parameters in the case of BQ Ser. The meaning of 
the symbols is the same as in Fig. 1. Top panel: the effect of 
setting T e ff = 5500 K; middle panel: setting incorrect gravity 
(logg = 3.0); bottom panel: incorrect microturbulence (v t = 1.5 
kms -1 ). 



5184 A , which is often used for radial velocity studies of late- 
type stars, was selected for this purpose. The cross-correlation 
function (CCF) was computed with the fxcor task. Having a rest- 
frame corrected spectrum, the EWs of the selected lines were de- 
rived with the task splot by fitting Voigt-profiles, allowing slight 
wavelength shifts. Those lines that showed shifts greater than 
+0.2 A from their laboratory position were omitted from the 
later analysis. 

For each line the iron abundances were determined with 
the SPECTRUM] spectral synthesis code by R.O. Gray. 



http://www.phys.appstate.edu/spectrum/spectrum.html 
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Table 2. Atmospheric parameters and iron abundances for the program stars 



Star 


Date 


S/N 




l°g£ 


logg + a 


v t 


A(Fe) 


(X 


(v e sin i) 








(K) 


(dex) 


(dex) 


(km s- 1 ) 


(dex) 


(dex) 


(km s" 1 ) 


Y Car 


Mav 30 
ividy j\j 


1 87 
1 o / 


6000 


9 

Z..VJ 


1 09 


4.1 


7 4S 


1Q 


u 




June 1 


270 


7000 


2.5 


2.48 


4.1 


7.43 


0.14 


10 


FY Car 








1 




4 9 


7 30 
/ .jy 


n ?4 






June 1 


139 


5750 


1.5 


- 


4.2 


7.24 


0.20 


20 


GZ Car 


IVldy J\J 


1 S8 






9 f)(S 


5 n 


7 4Q 
/ .^y 


0.22 


n 

u 




June 1 


125 


6000 


2.0 


1.98 


5.0 


7.35 


0.19 







ividy ju 




fiooo 

uuuu 


1 


? Of, 


3 8 


7 31 


n ? i 


A 




June 1 


168 


5750 


1.5 


1.59 


4.2 


7.29 


0.14 


8 


BKCen 


May 30 


201 


6000 


2.0 


2.08 


4.0 


7.65 


0.16 


10 




June 1 


151 


6000 


2.0 


2.17 


4.0 


7.48 


0.23 


4 


V1 048 f"Vn 


A/fav 30 

IVldy JU 


1 U 1 


69^0 

UZ JU 


z.u 




4. 9 


7 37 


n 93 


7 




June 1 


189 


6250 


2.0 


- 


4.2 


7.25 


0.22 


5 


VX Pup 


June 1 


274 


6500 


2.5 


2.51 


4.4 


7.29 


0.14 


10 


EWSct 


May 30 


150 


5750 


1.5 


1.57 


3.8 


7.48 


0.16 


10 




Mav 31 

IVldy .? 1 


1Q7 




1 <; 


n no 


3 8 


7 ^4 


1 8 


1 




June 1 


228 


6250 


2.0 


2.04 


3.6 


7.43 


0.20 


10 


V367 Set 


May 30 


151 


6000 


1.5 


1.33 


4.5 


7.54 


0.21 


15 




Mav 31 
ividy 3 i 


146 




? 

Z..V7 


Z..UZ 




7 

/ . jy 


n 9? 


i j 


V458 Set 


May 30 


218 


6250 


2.0 


1.98 


4.2 


7.55 


0.19 


16 




Mav 31 

IVldy J 1 


1 83 

1 J 


69^0 


Z.. J 


9 S4 


4 2 


7 (SI 


17 

U. 1 / 


1 j 


BQ Ser 


May 30 


133 


5750 


2.0 


1.81 


3.7 


7.34 


0.21 


17 




Mav 31 


1 30 




? 


1 07 

1.7 / 




7 30 
/ .jy 


n 97 


1 f, 

1 u 




June 1 


169 


6000 


2.0 


2.03 


3.7 


7.36 


0.25 


14 


U Tr A 


Mav 30 


300 






2.04 


4 8 


7 43 


14 






June 1 


270 


6000 


2.0 


2.33 


4.8 


7.25 


0.18 





AP Vel 


May 30 


126 


5750 


1.5 


0.00 


3.9 


7.51 


0.17 


18 




June 1 


154 


6000 


2.0 


2.01 


3.9 


7.35 


0.20 


10 


AX Vel 


May 30 


328 


6500 


2.5 


2.52 


5.1 


7.54 


0.17 


10 




June 1 


232 


6250 


2.0 


2.17 


5.1 


7.26 


0.19 


13 


GSC 8607-0608 


May 30 


169 


5750 


1.5 


1.50 


3.7 


7.38 


0.22 


7 


GSC 8691-1294 


May 31 


66 


6250 


2.5 




3.6 


7.48 


0.30 


15 




June 1 


274 


6000 


2.0 




3.3 


7.46 


0.15 


15 



SPECTRUM solves the radiative transfer equation in a plane- 
parallel model atmosphere assuming LTE. We used Kurucz 
model atmospheres with T e s = 5000-7000 K, logg = 0-2.5 dex 
as appropriate for Cepheids, and a solar abundance pattern. After 
selecting a particular atmosphere (specifying T e s and log g) and 
microturbulence v t , the LTE abundances A(Fe) = log(Np e /NH)+ 
12, where Np e and Nn is the number density of iron and hy- 
drogen, respectively, were calculated with the auxiliary code 
BLACKWEL for each input line. To minimize the number of 
lines that are affected by NLTE-effects, only weak lines with EW 
< 0.15 A were used. 

The atmospheric parameters were determined entirely from 
the observed EWs in the standard way (e.g. Fry & Carney [T997l l. 
The effective temperature is selected by requiring that the cal- 
culated abundances do not depend on the line excitation poten- 



tial. Similarly, the surface gravity is chosen when the calculated 
abundances are the same for both Fe I and Fe II lines (i.e. the 
abundances do not depend on the ionization state). Both T e ff and 
log g were further constrained by fitting the profiles of Ho- with 
model spectra calculated with ATLAS9/SYNTHE. Finally, the 
microturbulent velocity is set when the line abundances are in- 
dependent of the line EWs. 

Figure 1 illustrates the correct choice of these parameters in 
the case of BQ Ser. The final iron abundance (indicated by the 
horizontal line) was calculated as the average of all individual 
abundances including neutral and ionized Fe lines. The top, mid- 
dle, and bottom panels of Fig. 2 show the effect of an incorrect 
temperature, gravity, and v t , respectively. Decreasing either T e g 
from 6000 K to 5500 K (top panel) or the microturbulence from 
3.7 to 1.5 km s _1 (bottom panel), a strong dependence of the Fe 
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Fig. 3. The fitting of the Ha profile of BQ Ser. The observed 
spectrum is shown by continuous curve, while the model spec- 
trum is plotted as a dashed curve. The model shown here was 
computed assuming T e g = 6000 K and logg = 2.0. The fit was 
restricted to the wings of Ha, since the core cannot be modeled 
under the LTE assumption. 



abundance on the excitation potential is introduced. A change 
in logg from 2.0 to 3.0 (middle panel) causes a systematic dif- 
ference between the abundances calculated from Fe I and Fe II 
lines. 

One example of the Ha profile fitting is plotted in Fig. 3. 
Due to its temperature and pressure sensitivity, the Ha line is a 
good indicator of the photospheric temperature and gravity, even 
though the line core is more affected by NLTE effects. Thus, the 
fitting is restricted only to the line wings. Since the abundance 
data are only weakly sensitive to the gravity (see Fig. 2), the Ha 
fitting provided important constraints for this parameter. 

To test the abundance scale of BLACKWEL, as well as 
the accuracy of our EW measurements, the daytime solar spec- 
tra were also analyzed using the same line list. The Kurucz 
solar model atmosphere and its parameters (T e ff = 5777 K, 
log g = 4.4377, v t = 1 .5 km s _1 ) were used and the resulting line 
abundances averaged. The results were A (Fe) = 7.44 + 0.18, 
7.52 + 0.18, and 7.51 ± 0.17 for 30 May, 31 May, and 1 June, 
respectively, where the uncertainties were estimated as the stan- 
dard deviation of the line abundances around the mean. It is vis- 
ible that for the spectrum obtained on the first day (30 May) the 
calculated abundance is smaller than those from the two other 
days, but they agree well within lcr. Thus, we adopted their un- 
weighted average A (Fe) = 7.49 + 0.18 as the final solar abun- 
dance. This value is in remarkable agreement with the recent 
solar iron abundance of 7.45 + 0.05 (e.g. Asplund et al. 2005), 
and A (Fe) = 7.49 was adopted as the observed zero point for 
determining [Fe/H] from the FEROS spectra. 



4. Results 

The iron abundances, together with the corresponding atmo- 
spheric parameters, are collected in Table 2. The uncertainties 
of r e ff and \ogg are +250 K and ±0.25 dex, respectively, due 
the resolution of the applied grid of model atmospheres. The 
estimated uncertainty of the microturbulent velocities is +0.5 
kms -1 . The errors of the abundances are purely statistical, be- 
ing the rms standard deviation of the individual line abundances 
around the mean value (see the previous section). 



It is important to note that the inferred surface gravities are 
actually the effective gravities containing the acceleration term 
due to the pulsation 

GM dv mA 

where M is the stellar mass, R(t) the instantaneous radius of the 
photosphere, p the projection factor to convert the radial veloc- 
ities into the pulsational velocities (we used p = 1.38 from the 
period-dependent formula by Gieren et al. 1999), and v ra d the 
radial velocity at the moment of the observation. 

To get information on the true gravities of the program stars, 
we estimated the acceleration terms from the derivative of the 
radial velocity curves. The radial velocity data were collected 
from Stobie & Balona ( [19791 , Gorynya et al. ([19921 [19961 [19981 , 
Bersier et al. d!994| l, Antipin et al. dl999l l, and Petterson et 
al. ([2004). For each Cepheid, the radial velocities were fitted 
by a Fourier-polynomial including the /o, /i, /o + /i, /i - fo, 
and 2/o + f\ frequency components. The frequencies were com- 
puted from the periods listed in Table 1. Because these are de- 
termined from observations distributed on a long time baseline, 
they can be applied to derive the phases of the new observations 
unambiguously. Then, the analytic derivative of the Fourier- 
polynomial was computed at the phase of the given spectro- 
scopic observation to get the acceleration term in Eq. 1 . Table 2 
also contains the true gravities, i.e. the observed (effective) grav- 
ities corrected for the acceleration. 

The chemical abundances should not depend on time, so the 
data in Table 2 can also be used to test the stability and the pos- 
sible systematic errors of our abundance analysis. Generally, it 
is visible that the calculated abundances do agree within their 
errors. However, the spectra from day 3 (June 1) tend to give 
slightly lower abundances (but within lcr) than those from day 1 
(May 30). No such tendency can be seen in the atmospheric pa- 
rameters. We attribute this slight inconsistency to an instrumen- 
tal calibration problem, probably associated with the scattered 
light + sky background removal that may alter the line depths 
of continuum-normalized spectra. Nevertheless, since the differ- 
ences are not significant, we define the final iron abundances 
(Table 3) as the average of the values given in Table 2. Their 
final estimated uncertainties reflect their individual errors, plus 
the differences between the data belonging to different nights. In 
Table 3 we also list the new [Fe/H] values calculated from 

[Fe/H] = Aceph(Fe) - A (Fe) (2) 

where A (Fe) = 7.49 was applied (see the previous section). 
We note that 14 out of 17 double-mode Cepheids show slightly 
negative [Fe/H] values, though while considering their errorbars, 
the deviation from solar abundance is marginal. 

We checked our new iron abundances further by computing 
synthesized spectra with ATLAS9 using Kurucz models with so- 
lar abundance pattern scaled to the abundances in Table 2. Note 
that the solar iron abundance applied in the Kurucz models cor- 
responds to 7.63 on our scale, so the "abundance scale" param- 
eter was set to 10 [A(Fe) ~ 7 - 63 l, where A(Fe) is the iron abundance 
of the Cepheid. Moreover, we enabled the pulsation correction 
in VFIELD with the appropriate velocity value determined pre- 
viously. The final spectra are convolved with the FEROS instru- 
mental profile derived from fitting Gaussians to the spectral lamp 
emission lines. In Fig. 4 a synthesized spectrum is compared 
with an observed spectrum of EW Scuti. The agreement is satis- 
factory, although a few weak lines are not fitted well. The cause 
of this is that the log gf values in the Kurucz line list were not 
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Table 3. Final (averaged) abundances and uncertainties with [Fe/H] data from previous studies. For the latter data, uncertainties 
(where available) are given in parentheses. 



Star 


[Fe/H] 


cr 


[Fe/H] i 


[Fe/H] 2 


[Fe/H] 3 


[Fe/H] 4 




(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


YCar 


-0.05 


0.16 


-0.33 (0.35) 


- 


- 


- 


EY Car 


-0.17 


0.22 


- 


- 


— 


- 


GZCar 


-0.07 


0.20 


-0.36 (0.33) 


-0.08 


- 


- 


UZCen 


-0.19 


0.17 


-0.30 (0.32) 


-0.09 


- 


- 


BKCen 


0.07 


0.19 


-0.03 (0.30) 


— 


- 


- 


VXPup 


-0.20 


0.14 


-0.39 (0.14) 


-0.17 


-0.13 


-0.15 (0.12) 


EWSct 


-0.01 


0.18 


-0.08 (0.09) 


- 


+0.03 


- 


Vjo/ oct 


U.U/j 


n o 1 


+U.ZO (U.jU) 




— U.U1 




V458 Set 


0.09 


0.18 




- 




- 


BQ Ser 


-0.13 


0.24 


-0.36 (0.20) 


-0.08 


-0.05 




UTrA 


-0.15 


0.16 


-0.31 (0.34) 








AP Vel 


-0.06 


0.18 


-0.18(0.35) 


-0.045 






AX Vel 


-0.09 


0.18 


-0.43 (0.29) 


-0.12 




-0.15 (0.07) 


GSC 8607-0608 


-0.11 


0.22 










GSC 8691-1294 


-0.015 


0.22 










TU Cas 


-0.16 


0.24 


-0.43 (0.31) 


-0.05 


0.02 




VI 048 Cen 


-0.18 


0.22 











References to [Fe/H] values: 1: Andrievsky et al. J19941 ; 2: D'Cruz et al. ([2000 1; 3: Andrievsky et al. d20021 >; 4: Lemasle et al. (20071. 



1.05 




6150 6160 6170 6180 6190 6200 



Fig. 4. Comparison of the observed spectrum of EW Scuti (dotted curve) with a model spectrum computed by ATLAS9/S YNTHE 
(solid curve) with the parameters in Table 2. The solar abundance pattern is scaled to the calculated abundance of EW Set. Pulsation, 
as well as the instrumental broadening, has been taken into account (see text). 



adjusted to the values used by us in the abundance analysis (see 
above), which makes the strength of these lines different from 
those obtained with the SPECTRUM code using the same abun- 
dance. This test also illustrates that the iron line strengths and 



abundances computed by two independent codes, SPECTRUM 
and ATLAS9, are generally consistent. 

The new [Fe/H] data in Table 3 allow us to update the ob- 
served period ratio - the metallicity relation first pointed out by 
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o 
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[Fe/H] 

Fig. 5. The period ratio - metal abundance relation. The param- 
eters of the dotted line are given in Eq. 3. See text for the correc- 
tion term on the vertical axis. 

Andrievsky et al. ([1993 1994 ) in graphical form. From Tables 1 
and 3, the following empirical relation has been found: 

pjp = -0.0143 logPo -0.0265 [Fe/H] +0.7101 

+0.0025 +0.0044 +0.0014. {i) 

The period ratio, corrected for the weak dependence on the fun- 
damental period, is plotted against [Fe/H] in Fig. 5. 

4.1. The effect of rotation and pulsation on the derived 
abundances 

Rotation and/or pulsation alter the line profiles, thus, their ef- 
fect on the abundance analysis needs to be investigated. These 
macroscopic motions convolve the local line profile over the 
visible stellar surface, therefore the observed (integrated) line 
broadens and becomes shallower. As a first approximation, the 
total EW of the observed line does not change because of 
the conservation of the total absorbed flux. However, there are 
second-order effects, connected with e.g. the measurement of 
the EWs from the observed spectra with finite resolution and 
S/N, which may introduce some systematic uncertainties in the 
EWs of broadened lines, and these effects should be evaluated. 
We show below that in the case of our program stars neither ro- 
tation nor pulsation affects the results of the abundance analysis 
significantly. 

Cepheids are slow rotators, similar to other supergiants as 
expected from stellar evolution. The typical value of (v e sin z) is 
~ 10 km s _1 (De Medeiros et al. I2002.l t. There seems to be a 
correlation between the rotation and pulsation period: Cepheids 
with longer pulsation periods are slower rotators (Nardetto et al. 
2006). However, Takeda et al. ( fT997l l found (v e sin z) = for 
the short-period Cepheid SU Cas, which illustrates that even the 
shortest period (hence the smallest) Cepheids may be very slow 
rotators. Buchler and Szabo (2007) studied the effect of rotation 
on the period ratio of beat Cepheids from theoretical point of 
view, and find that the period ratio and its dependence on heavy 
element abundance is practically insensitive to the rotation (up 
to 20 km s" 1 ). 

To estimate the projected rotational velocities of our program 
stars, we determined the FWHM of a few weak, unblended lines 
(Fe I AA6056, 6096, 6157, 6165 A) and compared it with those of 
ATLAS9 model spectra (T etf = 5750 K, logg = 1.5, and v, = 4.2 
km s _1 ) rotated with 0, 5, 10, 15, and 20 km s _1 . We checked 
that variations in the atmospheric parameters up to 250 K, 0.5 



dex, and 0.5 km s _1 in the temperature, gravity, and microtur- 
bulence, respectively, do not alter the FWHMs of the rotation- 
ally broadened lines significantly. Thus, instead of computing 
detailed models for all program stars, we used only this sequence 
of theoretical model spectra as templates for estimating the rota- 
tional broadening. As for the earlier results mentioned above, all 
program stars were found to be slow rotators with (v e sin i) < 20 
km s _I (see Table 2). 

Since pulsation causes an additional broadening, using a 
non-pulsating template spectrum causes a systematic overesti- 
mation of the rotational velocities. Because of this, our estimated 
(v e sin z) values are actually upper limits to the true projected 
rotational velocities. To test the influence of pulsation on the 
model FWHMs, we added +20 km s _I pulsational velocity (typ- 
ical of our short-period Cepheids) to the model sequence and 
re-computed the FWHMs of the same lines. It was found that 
this pulsational velocity causes ~ 5 km s _1 overestimation of the 
(v e sin z), which is about the same as the uncertainty of the esti- 
mated rotational velocities (~ 3 km s~'). Note that the FWHM 
of the instrumental broadening profile of FEROS is ~ 0.15 A, 
which corresponds to ~ 7 km s _1 rotational velocity at 6000 A. 
Thus, (v e sin z) values lower than this level are not resolved well 
from our spectra. 

Finally, the effect of rotation and pulsation on the model EWs 
was found to be 1 % for v rot = 20 km s _1 and less than 0.1 % 
for Vp U i = 20 km s -1 . These numbers illustrate that the effect of 
rotation on the line EWs and the abundance analysis presented 
in the previous sections is negligible. 

Pulsation does not change the line EWs, but it does change 
the atmospheric parameters (T e ff, \ogg) along the pulsational cy- 
cle. However, the values reported in Table 2 show that the vari- 
ation in the derived abundances is typically less than the errors 
of the individual abundances, caused by mainly the uncertainties 
in the atomic data and the model atmospheres. Therefore, we 
conclude that, although for pulsating stars the determination of 
the atmospheric parameters is more complex than in the case of 
static stars, the analysis resulted in consistent iron abundances 
for different pulsational phases. Thus, the abundances for those 
stars that were observed only in a single phase are also probably 
not affected by the pulsational motion. 

5. Discussion 

In this paragraph we compare the spectroscopically determined 
physical parameters (r e ff, logg, [Fe/H]) of the program stars 
with results from previous observational studies and predictions 
from pulsation theory. In addition, we give an estimate of the 
color excess E(B - V) of the program stars from the spectro- 
scopic effective temperatures. 

5.1. Comparison with previous [Fe/H] determinations 

Among the stars listed in Table 3, four targets have [Fe/H] val- 
ues determined in this paper for the first time. However, for the 
majority of the program stars more than one value exists for iron 
abundance published previously, thus, it is important to compare 
our results with those of earlier studies. The results collected 
from the literature are also summarized in Table 3. In Fig. 6, the 
literature [Fe/H] values are plotted with respect to the ones de- 
termined in this paper. The uncertainties of our new [Fe/H] data 
are between 0.1 - 0.2 dex (Table 3). 

There seems to be a systematic difference of about lcr be- 
tween our data and those of Andrievsky et al. ( 119941 1 in the sense 
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that the latter show generally smaller iron abundance for the ma- 
jority of our beat Cepheid sample. However, the [Fe/H] zero 
point of Andrievsky et al. (1994), 7.51, agrees quite well with 
ours (7.49); thus, the cause of this discrepancy is probably the 
use of different input atomic data and Cepheid line EWs. Also, 
the spectral resolution of their data is lower than ours, so the 
corresponding errors are higher (~ 0.3 dex) than the uncertainty 
of the new data presented in this paper. On the other hand, the 
more recent [Fe/H] determinations by Andrievsky et al. (2002) 
agree much better with our data, although only 5 stars are com- 
mon in their more extended observational sample with ours. We 
believe that the agreement of the more recent (independent) data 
increases the probability that these are closer to the true [Fe/H] 
values of these stars than the earlier ones. Although Andrievsky 
et al. ( 120021 ) do not specify individual error bars for their [Fe/H] 
values, they note that the average uncertainties of their abun- 
dances are between 0.05 - 0.2 dex, which is comparable to the 
error bars of our data. 

The agreement is also good (within lcr) with the data by 
D'Cruz et al. (2000), who predict the [Fe/H] values from the- 
oretical pulsational models (and also point out the discrepancy 
between their results and those of Andrievsky et al. (1994). 
Unfortunately, they also do not assign errors to their theoreti- 
cal [Fe/H] values except a short note that "they should be quite 
accurate, provided the theoretical opacities and models are suffi- 
ciently accurate". Indeed, it is seen from Table 3 that they agree 
with our new spectroscopic abundances within the errors. 

Finally, only two stars are common in the sample of Lemasle 
et al. (120071 ) and ours. For both stars the results agree well with 
our ones and also with those of Andrievsky et al. (120021 ) . 

5.2. The metallicity of the 10/20 pulsator V1048 Cen 

The case of V 1048 Cen (=HD 304373) deserves a particular dis- 
cussion. Together with CO Aur, they are the only 10/20 pul- 
sators known in the Milky Way (Beltrame & Poretti 2002]). Its 
position in the Petersen diagram is very similar to that of other 
10/20 pulsators in the Large Magellanic Cloud. This similarity 
is striking, considering that the F/IO pulsators in the Milky Way 
and in the Magellanic Clouds describe well-separated sequences 
in the Petersen diagrams (see Fig. 1 in Poretti & Beltrame 2004). 
It could be argued that VI 048 Cen (and then CO Aur, too) 
could have a particular metallicity, different from those of other 
double-mode Cepheids. However, our data (see Table 3) demon- 
strate that the metallicity of VI 048 Cen is very similar to those of 
galactic F/IO pulsators. Therefore, it still remains unexplained 
why the two classes of double-mode pulsators behave differently 
in different galaxies. 

5.3. Period-gravity relation 

The linear pulsation theory predicts a simple relation between 
the period of the fundamental mode Pq and surface gravity, 
namely, 

logg = 2.62 - 1.21 log P (4) 

(Kovacs, 2000). This can be inferred from a simple pulsation 
equation and blackbody assumption. Using the spectroscopic 
gravities listed in Table 2, it is possible to test this simple the- 
oretical relation with observations. 

Since log g in Eq. 4 refers to the mean gravity, averaged over 
the pulsational cycle, the gravities corrected for the acceleration 
term (Table 2) were averaged. Only those stars that had more 



than one spectra were considered, although the average of two 
gravities obtained at random pulsational phases may be quite off 
from the proper mean value. Those stars were also omitted for 
which there were not enough velocity data to compute the cor- 
rection for acceleration (see Sect. 3). This resulted in 10 stars in 
the subsample. Their averaged (corrected) gravities are plotted 
in Fig. 7 as a function of logPo. (The error bars are estimated 
as +0.25 dex, although this may be optimistic due to the small 
number of averaged data.) The linear least-square fit to the data 
has the parameters 

logg = 2.76 - 1.20 logP ( „ 

±0.33 +0.54 , {D> 

while the theoretical relation is from Eq. 4. It is seen that the 
observational data support the theoretical relation nicely. The 
period dependence is fully recovered with only a slight zero- 
point difference. This gives further confidence to the reality of 
the calculated atmospheric parameters and abundaces in Table 2. 
A similar diagram is presented by Andrievsky et al. (2002)) for 
a larger sample of Cepheids, but they plotted the uncorrected 
gravities (containing the acceleration term) against log Pq. 

5.4. The period ratio - metallicity relation 

It is striking that, based on high-resolution spectroscopic ob- 
servations, the period ratios of the majority of known galactic 
double-mode Cepheids depend mostly on their metal abundance 
and (weakly) on the pulsational period. Since the simultaneous 
excitation of two pulsational modes (i.e. double-mode pulsa- 
tion), and particularly the period ratio itself, strongly depends 
on other stellar properties (mass, mean effective temperature, lu- 
minosity) as well, the dominance of the metallicity dependence 
suggests that the observed galactic beat Cepheids should have 
nearly the same mass and effective temperature. Indeed, the the- 
oretical masses of several of these stars inferred by D'Cruz et al. 
(2000) show a very low scatter of ~ 0.3 M around the mean 
value of 4.2 Af . 

We have also found an agreement between the observed pe- 
riod ratios in the Magellanic Clouds and the extrapolation of the 
galactic relation toward lower [Fe/H] values. This may mean that 
the slope of galactic relation is not strongly biased by the small 
number of galactic Cepheids analyzed here. A more thorough 
study of the metallicity distribution of Magellanic Cloud beat 
Cepheids will be the topic of a subsequent paper. 

5.5. Determining the color excess from spectroscopic 
temperatures 

The effective temperatures obtained in the manner described 
above can be used for determining the color excess of the tar- 
get Cepheids. The unreddened B - V color index corresponding 
to the given T e ff value was taken from Table 15.7 in Cox (2000) 
by interpolating the relevant data. With this lack of simultane- 
ous photometry, the instantaneous 'observed' B - V color index 
was determined from Berdnikov's (2006) extensive data set of 
Cepheids which contains multicolor photoelectric observations 
for 13 stars discussed in this paper (among several hundred other 
Cepheids). When predicting the color index for the moment of 
spectroscopic observations, the frequencies corresponding to the 
periods listed in Table 1 were used, and three harmonics of the 
fundamental frequency and the sum and the difference of the 
frequencies of the two excited modes were taken into account. 
Thus a pair of reddened and unreddened B — V values was gener- 
ated for each moment of spectral observation. Because the color 
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Table 4. The E(B - V) color excess of the program stars with 
available BV photometry 



Star 


E(B - V) 


YCar 


0.175 


EYCar 


0.135 


GZCar 


0.423 


UZCen 


-0.023 


BKCen 


0.163 


VXPup 


0.241 


EWSct 


1.089 


V367 Set 


1.273 


V458 Set 


0.636 


BQ Ser 


0.854 


UTrA 


-0.043 


AP Vel 


0.371 


AX Vel 


0.261 



X -0.1 
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Fig. 6. Comparison of the [Fe/H] data of this paper with values 
from the literature (see Table 3 for tabulated data). The meaning 
of the symbols is as follows: filled circles: D'Cruz et al. {2000); 
open circles: Andrievsky et al. (1994); crosses: Andrievsky et 
al. (2002); asterisks: Lemasle et al. ((2007). The reason for the 
discrepancy between our values and those by Andrievsky et al. 
( 119941 ) is described in the text. 

excess should not vary during a pulsational cycle, the color ex- 
cesses were averaged for the individual beat Cepheids, and these 
E(B - V) values are listed in Table 4. The comparison of the 
color excess data with those obtained previously by others and 
our own independently derived other values (from infrared pho- 
tometry) will be discussed in a subsequent paper. The accuracy 
of our color excess data can be judged from the physically unre- 
alistic, slightly negative E(B - V) values obtained for the bright- 
est (i.e. nearest) stars in this sample. Both UZ Cen and U TrA 
can be slightly reddened, therefore the error of the color excess 
derived here may be about 0.05-0.06 mag. 

6. Conclusions 

New, accurate, homogeneous iron abundances have been derived 
for a large number of galactic double-mode Cepheids (17 out of 
23) from unprecedented high-resolution, high S/N echelle spec- 
tra. For 4 of the program stars these are the first [Fe/H] values in 
the literature. The abundance analysis, based on Kurucz model 
atmospheres and LTE assumption, was computed in the canon- 
ical way. Besides [Fe/H] values, the physical parameters T e g, 
log g (corrected for acceleration of the atmosphere), microturbu- 
lent velocity, and upper limit for the projected rotational velocity 
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Fig. 7. Period - surface gravity relation from observations and 
theory. The dots are the averaged, corrected surface gravities of 
the program stars (see text). The dotted line is the relation pre- 
dicted by simple pulsation theory (Eq. 4), while the solid line is 
the result of a least-square fit (Eq. 5). The theoretical relation is 
nicely recovered from the observations. 

(v e sin i) were also determined from the spectra. For a few stars 
in common, our new abundances are in good agreement with the 
results of other recent abundance analyses in the literature. 

The observed period ratios of the program stars show strong 
dependence on the iron abundance and only weak dependence 
on the fundamental period. This may suggest that the scatter of 
the known Cepheids in the Petersen diagram (P{/Pq vs. logPo) 
is mainly due to differences in their metallicity. 
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